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Abstract 

Surface plasmon polaritons (SSP), moving along a smooth curved 
interface between two isotropic media with slowly varying dielectric 
permittivities and magnetic permeabilities and supporting SSP, are 
studied theoretically. Solutions of Maxwell equations are investigated 
within a small wavelength limit in the boundary layer of the wave- 
length order near the surface. An explicit asymptotic formula for an 
EM wave traveling along geodesies on the surface is obtained. An ex- 
ponential factor reflects a distinction between the planar and curved 
cases. A curvature-dependent correction term in the exponent demon- 
strates a strong dependence on transverse curvature and a weak de- 
pendence on longitudinal curvature. The closer the parameters to the 
resonance case, the more pronounced this tendency. It is found that 
the attenuation of the SPP in the case of lossy media may be reduced 
by changing the curvature. 

If the signs of the magnetic permeability of the medium on both 
sides of the interface are opposite, the surface magnetic plasmon po- 
lariton may propagate. Its short-wavelength asymptotics is found. 
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Fabrication of nanoscale devices, applications to near-field optics, mi- 
croscopy, and spectroscopy stimulated investigations of surface plasmon po- 
lariton propagation (see, for example [TJ [2j |3j H]). Progress in metamaterials 
attracts attention to magnetic surface plasmons [5]. The aim of the present 
paper is to derive an asymptotic formula, which describes surface-curvature 
effects of the SPP. Papers on the SPP on surfaces with spherical and cylin- 
drical geometries have appeared in the last century jS] , [7J [SJ EJ EE] ; now an 
interest in the dependence of properties of the SPP on the surface curvature 
has been renewed. The propagation of the SPP on bends is studied in [IT]. 
[12]. Numerical methods of finding the SPP on an arbitrary adiabatically 
varying surface, using a curvilinear coordinate system, were developed in 
[T3] and [13]. The numerical study of the influence of the curvature on the 
SPP characteristics was performed also in [T5] and [TB] . 

We give here an explicit asymptotic formula for the solution of the Maxwell 
equation near the interface between two media with smoothly varying dielec- 
tric and magnetic permittivities. We use a curvilinear system associated with 
geodesies, as was used in [T7J [T5]. Then we rewrite the Maxwell equation in 
matrix form by analogy with [TO] and apply the boundary layer method [20J . 
The distinction from the planar case is in a slowly changing amplitude fac- 
tor and a Berry phase, which provided a correction to the phase speed. In 
the case of lossy media, we obtain a curvature-dependent correction to the 
attenuation constant. The reduction of losses, as compared with the planar 
case, can be predicted with the formula. An explanation of this fact in the 
case of cylindrical bend has been done in [12J. 

A formula for a surface magnetic plasmon is found with the help of the 
symmetry properties of Maxwell'a equations. 

We solve the Maxwell equations for a monochromatic electric E and mag- 
netic H fields of a frequency u 

rot E = ikfiH., 
rot H = — ikeE, k — ou/c, 

c is the speed of light; the dielectric permittivity e and the magnetic perme- 
ability fi have a jump on the interface between two media and change slowly 
in each medium. We introduce a curvilinear coordinate system (s, a, n) as- 
sociated with an interface, where s is the length along the geodesic curve on 
the surface, the coordinate a parameterizes distinct geodesies, the coordinate 
n is measured along a normal at every point of the surface. We introduce a 
new variable vector \I/ = {£, H) 1 , where £ and Ti are covariant components of 
the electromagnetic field £ = g"E, H = gH, g is the metric tensor. In terms 
of parameters of the surface, namely, the longitudinal radius of curvature p, 
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the radius of curvature in transverse direction p t , the divergence of geodesies 
on the surface h, and the torsion r, the metric tensor is expressed as follows: 



( (l + ^) 2 + r 2 n 2 -hr(2n + ^ + ^) o\ 
-hr(2n + ^ + ^) h 2 ((l + f t ) 2 + r 2 n 2 ) 

V o o 1/ 



We rewrite the Maxwell equations in matrix form by analogy with [21j,[19| 



+ ir n d n * = -i r s d s - v a d a 
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The investigation of the exactly solvable problem of diffraction on a cylin- 
der shows that the plasmon wave is concentrated in the boundary layer of 
the wavelength order near the surface. We make an assumption n = O(r), 
where k is the wavelength in vacuum. We suppose that the unit of the length 
is the typical radius of curvature p and k >> 1. We decompose the matrix a 
in A; -1 and obtain 
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if = diag(/i, /i). We seek the solution of ([1]) in the form 

s oo 
ik f S(s',a)ds' ^ — , . , .. 

* = e o $(s,a,f), $ = ^AT J $ (j) . (3) 

Substituting (E} into the Maxwell equations and equating terms of equal 
orders of k, we obtain 

(p(°)+^ n )$(°) = /3f^ ), (4) 

-«(f s a s -a)$(°». (5) 

The systems of equations split into two systems in the upper and lower 
media. The condition of the continuity of the tangential components of the 
electrical and magnetic fields on the interface between the media yields a 
relation between the coefficients of the solutions. After simplifications we 
obtain 

<S>^ = C(s,a)cp, (6) 

iK + Q §_ 

he+ ' ' he-i 



K+i _ a j. 010 z>0; 

* = { r e+ h£ : v (7) 

.^,0,-^,0,1,0) , l<0, 

k± = \/ /3 2 — e±/i±, and the square root with a positive real part is chosen. 
The condition of continuity of the first component of #^°- ) yields the relation 
K + /e + = — K_/e_, which can be satisfied if E\ and Ei have opposite signs. 
This condition is equivalent to the requirement 

2 (a*+£- — H-e + )e + e- 

c c i 

where /3 coincides with the propagation constant of a plasmon on the plane 
interface between two homogeneous media. When the moduli of E\ and e<i 
coincide, a plasmon resonance occurs, f3 and k± tend to infinity, and the field 
may be concentrated in a subwavelength layer near the interface. 
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The function C(s,a) of the coordinates on the surface should be defined 
from the equation fl5j). We note that if e and fi are real- valued, then the 
operators on the left-hand side and right-hand side of (CQ) are Hermitian if 
the scalar product (f , g) is defined by the standard way 



OO g 

f,g)= / £#(0fc(Q<«; 

k=l 

— oo 



here, the superscript * stands for complex conjugation. The solvability con- 
dition of d^J) has the form 

J ($(rf s ^ (o) )+ J ($ (o) ,w (o) )+ 

($(o) ) p(i)$(o)) = . (9) 

We note that ($(°) , T S ^°A and (*(°) , f a $(°A are average fluxes of energy 

through areas of infinite height and unit width, which are orthogonal to s 
and a directions, respectively. Upon substitution of fl§D into this equation, 
we obtain an expression for C(s, a): 

C(s,a) = D(cp,t s cpY 1/2 ex P L / ^'^^ dsY (10) 

where D may depend only on a. After calculations we obtain 

/?£+/! 1 



^-wXh-TyiM-^- (i2 » 

Gathering formulas (j3J), (jSJ), ffTUj) . f fTTj) . f TT2|) . we obtain the principal term of 
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the asymptotics 



<°> =A exp | ik J (3 ds' J exp J % J (- + — ) ds 



a , °t\ j\ J exp (— £>0 

exp £ < 0, 
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fc W (;^->; ) V4 , o], (14) 

V W-e+ - fi+e-) (ei - e+) / 

° = T72' ( 15 ) 

2( £ + (-£_)(/i_£+ -//+£_) (//+£+ -//-£-) ) 7 

at= if ^ + -^- y /2 (-^ + ) i/2 v (16) 



2 \/J-S + — fl+S- J (—£+—£_ 

Here, the e± stand for £ + if £ > and for £_ if £ < 0. We has written (j^J) 
assuming for simplicity that e and are real- valued. In lossy media operators 
in ([[J are non- Hermit ian. The solvability condition of (jl]) differs from In 
all scalar products, the function <&W standing on the left must be replaced by 
a function <fr(°) + , which solves an equation adjoint to ([5]). Analogously, the 
function cp standing on the left in scalar products in (TTOT) must be replaced 
by <p + . To find <J>(°)+ we should take $^°- ) corresponding to /?*, £*, and /i*. 
Upon calculations, we obtain the same result (1T31) . (JSJ), (IT51) . (IT5|) . 

The formulas can be simplified in the case /i + = /i_ = 1; then formulas 
(USD read 

£+ + £- -l / 1f7 \ 

(T = , - = , 0- t = (7 . (17) 

v £ +(- £ -) 

The amplitude factor A varies adiabatically, because ft,, e, and fi may 
vary. The first exponential term does not feel the curvature, it is the same 
as in the case of homogeneous media. The second exponential term is a 
correction term, which is curvature dependent. It is real- valued in the case 
of nonabsorbing medium, and it represents the Berry phase. Both formulas 
([16]) and (TTT1) demonstrate a strong dependence on the transverse curvature 
radius near the plasmon resonance e + + £_ = 0. Moreover, the correction 
has a stronger singularity than the propagation constant f3 near the plasmon 
resonance. The influence of the longitudinal radius p tends to zero, when the 
conditions are close to the plasmon resonance and /x+ = /x_ > according 
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to f|T7|) . If p + and /i_ differ and are positive, the influence of p may be 
important but does not have a resonance character. 

A surprising result is the possibility of reducing losses by changing the 
curvature in the case where we deal with lossy media. The correction terms 
become complex- valued. They contribute to the attenuation of the SPP. 
The curvature radii have different signs for concave and convex surfaces. 
Therefore the correction terms may give not only attenuation but an increase 
in the wave as well. This effect is in agreement with numerical results on 
the circular bend [IT], [12]. An explanation of this phenomenon in the case 
of curved boundary the metal-dielectric is given in [12]. There are two loss 
mechanisms: the bending loss through the radiation into the dielectric and 
absorption loss in the metal. The bend increases the loss because of the 
radiation into the dielectric and decrease it because the bend pushes the 
wave out of the metal where the wave travels with attenuation. These effects 
can be comparable and can result in reduction of losses comparing with a 
planar surface [12J. 

Qualitatively new results can be obtained in metamaterials where mag- 
netic permeabilities may have opposite signs. Then the surface magnetic 
plasmon (SMP) may exist. To obtain a formula for the surface magnetic 
plasmon, we should use the symmetry of Maxwell equations and interchange 
e and p, E and — H in the result ffT0]) - (fT6]) . We use the notation (3 m , a m and 
Otm for SMP instead of the notation (3, a, and a t for SSP. The propagation 
constant (3 m reads 

o2 Q-£+ ~ p+E-)p+P- 

Pm ~ o _ „2 • 

The dependence of the exponent in the correction term on p and p t is con- 
trolled by the parameters 



2 (p l+ (-p J J)(e-p J+ - e+n-) (e+p + - e-pJ) ) 1/2 
1 / s + p + — £_ / u_\ 1 ^ 2 (— / u_ / u + ) 1 / 2 



2 \e_p + — £ + /i_y (— p + — //_ 



If e_ and p- are negative, both types of plasmons may exist. The special 
cases where ps + = p + £- or p + e + = p^e^ correspond to the case of the 
coincidence of propagation constants (jHJ) and (ITS]) . In this case, both asymp- 
totic formulas are not applicable, plasmons do not propagate independently, 
and new formulas describing the process of plasmon interactions should be 
obtained. If parameters are close to degeneration but the formulas are still 
applicable, the correction terms may be very significant. 
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